The Luther-Emery liquid: Spin gap and anomalous flux period 
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We study the dependence of the ground state energy on an applied Aharonov-Bohm flux $ for the 
Luttinger model with large momentum scattering. Employing the method of finite size bosonization, 
we show that for systems with a spin gap but with gapless charge degrees of freedom, the ground 
state energy has an exact period of hc/2e, i. e. half a, flux quantum, in the limit of large system 
size L. Finite size corrections are found to vanish exponentially in L. This behavior is contrasted 
to that of the spin gapless case, for both even and odd particle number. Generalizations to finite 
temperature are also discussed. 
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I. INTRODUCTION 

Models of interacting electrons in one spatial dimen- 
sion are very valuable for the understanding of strongly 
correlated systems. This is because there exist theoretical 
methods enabling us to determine their physical proper- 
ties reliably. Indeed, by combining perturbative renor- 
malization group 1 , bosonization, and Bethe ansatz tech- 
niques, a wealth of interesting phases in one dimension 
has been discovered. 

While some properties of these phases are unique in one 
dimension, others have their higher dimensional analogs. 
For example the independent gapless spin- and charge- 
excitations and the vanishing quasiparticle weight of the 
Luttinger liqui$> are unique in ID. However, the fact that 
it has a finite charge compressibility and Drude weight is 
analogous to a normal metal in higher dimensions. As an- 
other example, like systems in higher dimensions, a Mott 
insulating state is realized at half filling for repulsive in- 
teractions. However, the fact that antiferromagnetic long 
range order is absent and that spin 1/2 excitation exists 
in the half-filled Mott state are special features of ID. 

Furthermore, in one dimension there exists a phase, 
the Luther- Emery liquid^, which exhibits a spin gap 
and no charge gap. In addition, as in the Luttinger 
liquid, the DC electric conductivity is infinite. The 
above characteristics suggest that the Luther-Emery 
liquid is a ID analog of a superconductor. However, 
up until very recently an important question remained 
unanswered: "Do electrons pair in the Luther-Emery 
liquid ?" The best way to answer that question is to 
determine whether the magnetic flux period is hc/e or 
/ic/2e4' 5 However, since the spin and charge degrees of 
freedom are manifestly separate in the effective theory 
describing the Luther-Emery liquid, and the vector 
potential enters only in the charge action, it is difficult 
to see why the flux period for a Luttinger liquid and a 
Luther- Emery liquid should be different. 

In a recent paper we addressed these issues in the one- 
dimensional t-J-J'-model in the limit of vanishing ex- 
change couplings— Fortunately, both a spin gapless phase 



as well as a spin gapful phase appear in this limit— In 
Ref. we have demonstrated that while the flux period 
is $o = hc/e in the former—, it indeed becomes hc/2e in 
the latter. In particular we have shown that as a function 
of the Aharonov-Bohm flux, the ground state energy of 
a spin gapped ring is periodic with period hc/2e. Due 
to one-dimensionality the energy barrier between adja- 
cent minima is proportional to the inverse circumference 
L of the ring. For definiteness, we therefore define the 
function 



£($)= lim L (£„(*)- #o(0)) 



(i) 



where Eo($>) is the ground state energy of the system as 
function of flux. 

Despite the above progress, the question "do all 
Luther-Emery liquids exhibit an hc/2e flux period, and 
hence electron pairing ?" remains to be answered. In this 
paper, we show that the answer to the above question is 
indeed affirmative. Technically we start from the Lut- 
tinger Hamiltonian with the g\ channel scattering— We 
bosonize this model using the constructive formalismSiifl 
which provides rigorous operator identities on the 
Hilbert space of the finite size system. We show that 
due to a set of constraints on the total charge/spin 
number/current operator a 2 ' 11 , the state of the spin sec- 
tor impacts the charge sector through a twisting of the 
boundary condition. As a result, when the spin sector 
is gapped by the large momentum transfer two body 
scattering, the charge channel flux period becomes hc/2e. 

In the literature, the fact that there exist constraints 
on the total charge/current operators in bosonization has 
been employed by Loss^ for spinless fermion systems 
to study particle number parity effects. Regarding 
spinful fermions, Ref. used a method similar to ours 
to determine the flux period for the Hubbard model. 
However, the author concluded that the flux period is 
always hc/2e regardless of whether a spin gap exists, 
which we believe to be in error. Furthermore, a common 
reasoning encountered in the literature is to attribute 
the hc/2e flux period to the dominance of singlet 
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superconducting (SS) correlations at long distance and 
low energy, rather than to the appearance of a spin 
gap. It has, however, been noted that states with 
dominant charge density wave (CDW) correlations may 
a lso fea ture this anomalous flux period (see, e. g., Refs. 

Here we argue that this is just the case when 
there is a spin gap. In this case, it is natural to interpret 
the state as being formed by Cooper pairs. The degree 
of coherence of these pairs will determine if the state is 
more appropriately thought of as CDW-like or SS-like on 
not too large length scales. In this picture, one naturally 
expects the flux period to be one half of a flux quantum. 
In the following, we will show that regardless of the 
correlation functions in the charge sector, the existence 
of a spin gap alone indeed causes the hc/2e flux period 
in systems with even particle number. 

The structure of this paper is as follows: In section 
ITT1 we present the Luttinger model with large momen- 
tum scattering and state the selection rules between 
charge and current quantum numbers that characterize 
its Hilbert space. In section ITTT1 we briefly review the for- 
malism of constructive bosonization and introduce some 
notation. In section IIVI we complete the proof that the 
flux period will be hc/2e in the presence of a spin gap, 
and contrast this behavior with that expected in the spin 
gapless case for even and odd particle number. We will 
also comment on finite temperature effects here. Our 
conclusions are summarized in section ]V\ Appendix lAl 
discusses the finite size refermionization of the spin part 
of the Hamiltonian, supplementing our line of arguments 
given in the bulk of this paper. Appendix iBlis devoted to 
the use of conjugate phase variables in the construction 
of Klein factors. 



II. THE MODEL AND THE SELECTION RULES 

The Tomonaga-Luttingcr Hamiltonian describes a gas 
consisting of right and left moving chiral fermions, each 
suffering small-momentum transfer scattering in a one- 
dimensional system of size L: 

Htl = Hq + Hi + H4 (2) 

H = ( rVpk ~ V) : C lks C rks ■ (3) 

r,k,s 

H 2 = j- Y (92\\8s, s > +92±5s,-s') p+. s (q)p-,s'{~q) (4) 

q,8, s f 

Ha = — Y (04||£>,s' +94x8s,-s>) ■ Pr,s(q)p r ,s'(-q) ■ 

r,q,s,s' 

(5) 

Here k = 2nn/L denotes the allowed momenta under pe- 
riodic boundary condition, the fermion operator c r k s an- 
nihilates a right (r — +) or left ('" = —) moving fermion 
with momentum k and spin s (see Fig. p = vftt/L, 




FIG. 1: Right- and left-moving branches of the Luttinger 
model. The crosses denote the allowed momenta k = 2irn/L 
for periodic boundary conditions. The dark shaded region 
represents the occupied momentum states in the "vacuum" 
N r . 3 — 0. The vacuum chemical potential lies between the 
last occupied and first unoccupied states as indicated by the 
horizontal line. The light shaded region corresponds to a dif- 
ferent filling. The Fermi momentum corresponding to the 
latter is given by kf = 2nN riS /L. 

and 

:O:=O-(O) . (6) 

In the above, < ... >o denotes the expectation value 
taken in the vacuum state defined as the ground state 
of Eq. ©. The density operators appearing in Eqs. (J2J- 
|JS| are defined as 

PrAq) = : c l,k+q,s C r,k, S ■ (7) 

k 

The q = -component of these operators, 

N r , s = p riS (0) (8) 

measures the extra number of (r, s) type fermions added 
on top of the vacuum. All four integers N rjS are conserved 
by Htl- These quantum numbers play an important 
role in the rest of the paper. Their importance in the 
bosonization procedure has been stressed by Heidenreich 
et ali& and HaldaneiS. 

Out of the four operators N r>s we can form the follow- 
ing linearly independent number and current operators: 

N P = J2 r>s N r , s , V . r.Y,.. 

N* = Er,s * N r,s , J* = Er, S ™ N r , s , (9) 

where the indices p and a stand for charge and spin 
respectively. It will be important in the following to 
note that in any one-band model with single particle 
states symmetrically occupied between kf = 2irN ryS /L 
and — kf, the total particle number is actually given by 

N = 2 + Y N r , s = 2 + N p (10) 

r,s 
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The reason for this is that the states at k — 0, which con- 
sist of 4 degenerate states in the Luttinger model rather 
than 2, have not been included in the definition of the 
N ryS (see Fig. [TJ. 

There are important relations between the integer 
quantum numbers defined in Eq. 0. For example 
N,N p a , Jp.a are either all odd or all even. In addition, 
the average of N p and J p has the same even-odd parity 
as the average of N a and J a , while they both have oppo- 
site even-odd parity as the average of N and J p . These 
constraints are summarized by the following "selection 
rules"ii: 

(-1)* - (-1)*' = {-l) J > = (-1)*" - {~l) J ° (Ha) 

- (-l)( Ar + J p)/ 2 = (_1)( W P + - 7 P)/ 2 = (_l)Wr+-M/2 

(lib) 

which follow from the definitions Eq. and the fact 
that the N r s are integer. For most of the paper, we 
shall primarily concentrate on the case where N is even. 
While selection rule Eq. I|lla|) then requires the same of 
all the other quantum numbers, it is the selection rule 
Eq. Ijllbjl that imposes a coupling between the spin and 
charge quantum numbers which ultimately determines 
the value of the flux period. 

The Tomonaga-Luttinger Hamiltonian Eq. (J2J is ex- 
actly solvable i2ii2ii& The solution describes a system with 
gapless spin and charge excitations. A spin gap may be 
opened by the addition of the following large-momentum 
transfer scattering term£: 

H = H TL + H 1 (12) 

X 1 C l,k'+q,s' C +,k,s ■■ C- tk - q , s C-,k',s> ■ (13) 

When the number of particles is incommensurate with 
the number of lattice sites, Eq. I|12|) is the generic Hamil- 
tonian including the most relevant two-body scattering 
terms. The inclusion of H\ destroys the exact solubility 
of the model, and at the same time it destroys the conser- 
vation of J u . However since Hi changes J a in multiples 
of 4, the parity (— 1) J ^/ 2 remains conserved. As a result 
the selection rules Eq. (|llaf> . Eq. ijllbfl remain valid even 
in the presence of Hi . 

III. BOSONIZATION 

Under suitable choices of parameters, Eq. (fT2*)l can 
describe a translationally invariant system of spin-1/2 
fcrmions with a spin gap but no charge gap, i.e., a Luther- 
Emery liquid. In the rest of the paper we study the de- 
pendence of the ground state energy of such a model as a 
function of an applied Aharonov-Bohm flux. Technically 



we employ the constructive bosonization method&iil ex- 
tensively reviewed in Refs. and In the following 
we shall just summarize the main bosonization rules. 

Due to the following commutation relation between the 
density operators 

kL 

[pr,s(-rk),p r i s >(r'k')} = —S rr >6 ss >6k,k> ■ (14) 

Ztt 

we define boson creation operators for each momentum 
q ^ and each spin s: 



[b s (q),bl,(q')} = 5 s , s ,S ml , [b s (q), b s ,(q')} =0,(15) 

where Q(x) is the Heaviside step function. The bosoniza- 
tion of the local fermion operators 

W,.(aO = -^fY.^ ^) ( 16 ) 

then proceeds by means of the introduction of a non- 
Hermitian bosonic field 




in terms of which the fermion creation operators can be 
written&ifiiil as: 

A rs (x) = -^^se^We'^We'^, (18) 

where the factor 

AriS =e i ^'' N -^ +s ^'' N ^'), (19) 

is introduced to ensure the proper anticommutation rela- 
tions between the fermion operators H18[) carrying differ- 
ent r and s. It commutes with all the spatially dependent 
fields in Eq. (JTSJ. A positive infinitesimal a was intro- 
duced in Eq. 1|17|) to ensure the convergence of commu- 
tators between operators. The operator Tp r s is conjugate 
to N ns , 

[7p riS ,N r>s ]=i. (20) 

Note that the validity of Eq. I|2U|I formally requires N riS 
to have a continuous spectrum (see Appendix lB|l . This is 
clearly not the case in the physical Hilbert space H p h ys 
we have been working in so far. We find it convenient, 
however, to introduce a larger Hilbert space 7i, where the 
N r , s operators have a continuous spectrum. This con- 
struction is analogous to the embedding of a discrete lat- 
tice into a continuous space, and is reviewed in Appendix 
IbI To ensure that the Hamiltonian, as well as physical 
observables, do not lead out of Ti. p h ys , the operators Tp r s 
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may only enter through integer powers of the unitary op- 
erators exp(z^ r s ) , which raise the N r>s by 1. We shall 
have occasion though, e. g. in Appendix 1X1 to work in a 
larger subspace of TL defined below. Formally it is most 
convenient to define operators that are valid everywhere 
in TL. The anticommuting operators A r ^ s exp(itp r s ) are 
also known as Klein factors in the literature. 

It is customary to further define local Hermitian fields 
each associated with the spin (a) or charge (p) degrees 
of freedom, 

(f)p.a(x) = - r (<Pr,l(x) ± fr,l(x) + h.C.) + p , CT 

r 

^ = iS r ^ ± ^' (21) 

r 

as well as their "dual" fields, 



4 

1 

4 ^ 



(22) 



Written in terms of the spin and charge boson operators 



1 



b p ^(q) = ^=(b^q)±b l (q)), 
the above local fields read (y = p, a) 



(23) 



— 7T N V X i ^ / 7f 

"M 35 ) = ?V - 77— F 2 Z^ s S n W 



2 L 2 ^ W V £|<?| 
x e-^-l^l"/ 2 (ftt + (24) 



x e-^-l'W 2 (6t ( 9 ) - 6 „(_ g )) . (25) 

From Eqs. (|24() and (|25|l it is evident that N v and J„ are 
the winding numbers of <j) v and 9 V respectively, and <f> v 
and B v are the spatial averages of <j) v and It is simple 
to check that <$> v and are the conjugate operators of 
J v and N v respectively, i.e., 



$ v ,j v ] = i, [e v ,N v ] = 



(26) 



Note that although the commutation relations (|26[) are 
analogous to Eq. J20J, the operators e 1 * 1 - and e l0u lead 
out of the physical subspace. This is so since within 
7~tphys-i the quantum numbers N v , J v cannot be raised or 
lowered by 1 independently, but are subjected to the se- 
lection rules Eq. iJTTJ. Within this space, only powers of 
e 4l< /v anc j e 4«e„ are allowed. 1 - However, within the larger 
space TL introduced above, the operators e % ^» and e %e " 
are nonetheless well defined objects. It is convenient to 
introduce a space of "fractional" excitations, TLf rac , gen- 
erated by acting on TLphys with all possible combinations 



of e % ®» , e %6v . Within TLfrac, the quantum numbers N„, J v 
are independent integers. We must bear in mind, though, 
that all physically acceptable states live in TLphys- 

The inclusion of the zero modes <j) v and B v in Eqs. 
(|24|l and 1251) ensures the proper commutation relations 
of these fields when the system size L is finite 

y>„(x),e u >(x')] =i^8 u ,u'Sgn(x-x'). (27) 

Eq. (|2*7|l suggests that the conjugate operator of 4> v {x) is 
proportional to d x 6 ll (x), i.e, 



U u (x) = d x 9 v {x) 

TT 

[<l> v (x), H„i (x')] = i 5 u y S(x — x'). 



(28) 



Similarly the conjugate operator of 8 v (x) is proportional 
to d x (j} v {x), i.e., 



v(x) 



■d x (j> v (x) 



[8 v (x),v'(x')]=i8 v yS(x-x'). 



(29) 
(30) 



The physical spin or charge density is given by Eq. Il29|l , 
and in the absence of an applied flux, the physical (spin 
or charge) current density is given by 



j u (x) = K v v v TL v {x) = K v v v d x 8 v (x), (31) 

which follows from the bosonized Hamiltonian given be- 
low. In the above expressions 



9Av\ 2 (92v^ 2 



v F 



K v = 

9lp,a = 
9ip.(J ~ 



/ttvf - 


" 92v - 


f 9Av 


ttvf - 


f 92v ~ 


f 9iv 


92\\ ± 


92± - 


9i\\ 




2 




9q\ ± 


94-L 




2 







(32) 



Since H p is just the density of right moving fermions mi- 
nus the density of left moving fermions, it is appropriate 
to interpret the coefficient K p v p = v F in Eq. JsHJ as the 
renormalized Fermi velocity of the system. 

In terms of <j> v {x) and 8 u (x) the bosonization identity 
Eq. JTBJ) reads 

iPl s (x) = — A rs : e i ( s p( x )+ r <t>p(x)+s( °( x )+ r <t>-( x ))) ; 

(33) 

Here, : () : denotes boson normal ordering: all powers of 
the fields ipl are to be moved to the left of powers of 
the fields f r ,si whereas positive powers of the operator 
exp(i<p rs ) are to appear on the very right, and negative 
powers of the same operator are to appear on the very 
left of the expression. 
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By means of Eq. I|33[) the selection rules Eq. be- 
come equivalent to the requirement 

ipr,s(x) — ■>pr,s( x + L) Vr, s. (34) 

This clearly illustrates the topological origin of these 
rules. 

We may now write the Hamiltonian Eq. i|12|) entirely 
in terms of the bosonic fields introduced above. The 
Tomonaga-Luttinger part of the Hamiltonian, including 
the large momentum scattering term with parallel spin, 
takes the following quadratic form: 

Htl + -Hi, 1 1 

= E v » \ E M WAq) + ^ {K/Ku + Jl K v ) 

(35) 

= E 7 / dx ■\Ku{d x 9 v {x)) 2 + ^L(d x <j> v {x)) 2 Y 

(36) 

where the operators b„(q) are related to those in Eq. ll'-'il) 
by a Bogoliubov transformation. The large momentum- 
transfer scattering with antiparallel spin term becomes^: 

#i,± = -^ [dx :coB(4^ ff (ar)): (37) 

Note that the coefficient of 4 in the argument of the cosine 
assures that the operator does not go out of the physical 
subspace, as explained above. 

The weak coupling renormalization group flow of the 
system Eq. <jT2|) is well knowriSS: For K a < 1, the op- 
erator i?i,j_ is relevant and a spin gap will be opened. 
This is the case we will focus on in the following. For 
spin SU{2) invariant systems <?ji = gtj_ = gi. In that 
case K a < 1 requires g\ < 0, as discussed by Luther and 
Emery. 3 



IV. THE FLUX PERIOD 
A. The spin gapped case 

By virtue of Eqs. (J3SJ) and jSZJl, the model Eq. JEJ) 
takes the form 

H = H TL + H 1 = H p + H a (38) 

where H p and H a act exclusively on charge- and spin- 
degrees of freedom respectively. The eigenstates are thus 
direct products of charge states and spin states 



and the ground state energy is the sum of spin and charge 
energies 

E Q = E C + E°. (40) 

When Hi causes a spin gap to open up, the spin sector of 
the model Eq. (|12fl is described by a gapped sine-Gordon 
field theory. 

In the following we shall focus on the N a — sector, 
which is where the gapped spin ground state lies. In 
this sector Hi j_ is relevant, and one may interpret the 
cosine term in Eq. I|37|l as a steep potential experienced 
by tp a (x). In the limit of infinite system size where true 
symmetry breaking is possible, one may think of 4> a (x) 
as being locked to one of the minima of the cosine po- 
tential. When this happens </> CT , the spatial average of 
<j) v {x), will take a c- number value equal to the respec- 
tive minimum value of cf) a . At first, let us neglect the 
selection rules Eq. 1111) . That is, we start by looking at 
the problem in the space 7if rac introduced in Sect ion ITTT1 
where in particular J a is an independent integer valued 
quantum number. Then we may regard the conjugate 
variable <j) (J as an angular variable with period 2ir. This 
notion becomes precise if we identify (j) a with its "lattice 
version" discussed in Appendix [51 which we shall do for 
the present purposei^i Within [0, 2ir) there are four in- 
equivalent minima of the cosine term in Eq. I|37|l . and 
the corresponding ground states in the spin sector can 
be labeled as 

|0} , lvr/2), |tt), |3tt/2), (41) 

where 

^\<f>)=<f>\4>)- (42) 

As discussed earlier, the operator 

fj = (~1) J ^ 2 = exp (-mtJ ct /2) (43) 

commutes with H Gl hence its eigenvalues can be used to 
classify the spin ground states. Unfortunately the states 
given in Eq. (|41Jl are not eigenstates of t). Following 
Appendix iBl it is easy to show that 

^fj\4>} = (0+ 1) m> (44) 

where the eigenvalue on the right hand side is to be un- 
derstood modulo 2tt. We may hence choose the global 
phases in Eq. (|41|l such that 

Iztt/2) = rf |0) . (45) 

It is thus easy to form linear combinations 

3 

l*7>=X>~**n o > (46) 

such that 



|c) ® \s) , 



(39) 



v\v) =v\v) 



(47) 
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We are now in a position to enforce the selection rules 
(|llfl . Given N a = 0, the selection rule Eq. (|lla|) requires 
Jcr to be even. As a result only rj = ±1 are allowed. We 
label these two states by 

I-)- (48) 

Thus actually, the ground state is only two-fold degener- 
ate. This degeneracy becomes further lifted in the case 
of a finite system size L, to be discussed next. 

For finite L, the notion that the field <fr a is locked to 
a classical value is no longer valid. In fact for finite L, 
even <\> a is subjected to quantum fluctuations. This is 
explicit in Eq. I|35[) . where the variable conjugate to 4> a , 
namely J CT , enters the Hamiltonian when L is finite. Thus 
the spin ground state can no longer be thought of as one 
of the "locked" spin states given by Eq. I|41|l . On the 
other hand, since rj remains a good quantum number, 
the states in Eq. I(48|l are still well defined as the re- 
spective ground states in the rj = ± sectors of the spin 
Hilbert space. We note that the spin states in Eq. 148(1 
thus defined are not strictly degenerate for finite L. ft is 
important to observe, however, that the difference in en- 
ergy between these two states vanishes exponentially in 
the system size L. One way to see this is the well known 
fact that the gapped sine-Gordon field theory is the low 
energy effective theory of a dimerized spin- 1/2 chainiS 2 . 
Here, the |±) are respectively the symmetric and anti- 
symmetric combination of the two dimer patterns. Since 
the two dimer patterns differ by a macroscopic number 
of degrees of freedom, the tunnel splitting between these 
two states should vanish exponentially with the system 
size. A slightly more direct way to see the above is of- 
fered by the well known mapping between the gapped 
sine-Gordon theory and the massive Thirring model, 3,23 
We will elaborate on this point in Appendix [S] The 
advantage of this method is that at the special Luther- 
Emery point, it allows us to study the effect of a finite 
temperature. 

For the purpose of this paper we may ignore the above 
exponentially small energy difference between the states 
Eq. I|48|l . This is because such a tiny difference will drop 
out of in the limit taken in Eq. In this sense we 

may still speak of a degeneracy in the spin sector of the 
model, and regard the spin contribution Eq in Eq. l|4l7|) 
as essentially independent of rj in the spin gapped case. 

Naively the spin degeneracy discussed above seems to 
suggest that the ground state of the full Hamiltonian 
Eq. (|38[1 is degenerate. However this is not so, and the 
reason for this is the selection rule Eq. ijllbfl . To demon- 
strate that let us assume the total particle number to be 
TV = Am + 2, whereas N a — 0. According to the selection 
rule Eq. Ijllb|) the spin states |±) may not be combined 
with the same charge state. The spin state |— ) may only 
be combined with a charge state whose current quantum 
number J p is an odd multiple of 2 and hence non-zero. 
The presence of a non-zero current will cost an energy of 
order v* F /L as is evident from Eq. (|3"5|) . The state |+), 
on the other hand, may be combined with a charge state 



of zero current, which minimizes the charge energy. As 
a result there is an energy splitting ~ l/L between the 
lowest energy state in the rj = + and ij = — sectors. We 
note that an analogous result was discussed by Haldanc 11 
for the case of a vanishing spin gap and a finite charge 
gap at commensurate band fillings. In contrast, here we 
are interested in the effect of an applied Aharonov-Bohm 
(AB) flux, which is of interest only when the charge sec- 
tor is gapless. 

The coupling to a vector potential A(x) is determined 
by gauge invariance and can be worked out from the mini- 
mal coupling requirement. We only consider the constant 
vector potential A(x) — Q/L corresponding to an AB 
flux. The correct coupling to <E> then follows from the 
formal replacement^ 

rl>l t ,{x) —> e~' l ^^ x ^ rs (x) (49) 

in the Hamiltonian, where a charge — e is assumed. Here, 
the boundary conditions of the field ipl s (x) remain the 
same, while the right hand side of Eq. (|49|l will in general 
satisfy different boundary conditions. By Eq. (|33|) . this 
is equivalent to the following replacement in the Hamil- 
tonian Eq. I|3(j|) and the current Eq. I|31l) : 

e p {x) — > e p (x)-^^ x (50) 

or, by Eq. I|25|l . simply 

J p ^J p +4§-. (51) 

Note that we did not attempt to introduce the gauge 
flux prior to bosonization. This is due to the fact that 
the fermionic field theory @ suffers from the well known 
chiral anomaly. The latter renders the global current of 
the model ambiguous in the presence of a general AB flux, 
unless gauge invariance is manifestly enforced. Through 
the "Lenz rule" / = -cdE (§) /<93> (Kohn, Ref. H see 
also Eq. (G2J), this ambiguity also enters the ground state 
energy dependence on flux. To deal with this problem 
will in any case require the use of gauge invariance, and 
this is most conveniently achieved in the final, bosonic 
language. It would be interesting to obtain Eqs. JSUJ, 
(|51|l via a more "microscopic" route, i. e. via bosoniza- 
tion of a microscopic (lattice) Hamiltonian with flux; this 
is subject to current investigations. We stress again, how- 
ever, that Eq. I|5U|) is uniquely determined by the minimal 
coupling principle. 

From Eqs. i|35[l . <|51l> the energy versus flux function 
£(<!>) in Eq. Q is given by 

7T / $ \ 2 

£ ($) = min - v p K p [J p + 4 — (52) 

J p =. ..-2,0,2... 4 1 ' \ H <f>o/ 

Here, all multiples of 2 are allowed values for J p by se- 
lection rule Eq. (jllafl . This leads to the various branches 
shown in Fig. [21 The alternating labels of rj — ± reflect 
the fact that the spin state has to be adjusted according 
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to selection rule Eq. (|llbjl whenever J p is changed by 2. 
In the presence of a spin gap, however, this does not af- 
fect the energy in the limit of Eq. 0J, as discussed above. 
As a consequence, Eq. (15211 has an exact period equal to 
half a flux quantum, shown by the lower envelope in Fig. 
Eli). We note that these findings are in complete agree- 
ment with those obtained in Ref. |y for the t-J-J' model. 
The amplitude of the ground state energy modulations 
is apparently given by 



AE 



AL 



-K v 



AL 



(53) 



where v* F is the renormalized Fermi velocity introduced 
above. The corresponding modulations of the charge cur- 
rent for a given quantum number J p are given by 



h = 



-e-jr J dxj p {x) 
j^ v f J dx (d x 6u(x) 



2tt $ 
L $q 



L 



$ 



(54) 



The current is thus diamagnetic for — *E»o /2 < $ < $o/2 
and is given by a sawtooth curve in general which one ob- 
tains by taking the derivative of the envelope in Fig. |2k>). 



The amplitude of the current is given by AI 



r/L, 



which is the same as that of spinless particles 12 , although 
the flux period is halved. Note that this observation is 
consistent with the notion that the charge of the carriers 
is effectively doubled. 



B. The spin gap less case 

To establish the fact that the hc/2e flux period is due 
to the presence of a spin gap, it is prudent to demonstrate 
the change of flux period when the spin gap collapses. 
First, let us assume N = 4m+2 as before. Our discussion 
from the preceding sections generalizes most easily to the 
case of a vanishing spin gap, if we also assume isotropy in 
the spin sector: In this case, SU(2) invariance requires the 
parameter K a to be unity at the Luttinger liquid fixed 
point. We will comment on the general non- isotropic case 
below. 

For gapless spins, the operator in Eq. (|38[) is ir- 

relevant, and we may expect to get qualitatively correct 
results by omitting it. With this simplification, the spin 
sector becomes analogous to the charge sector, and in 
particular J a can be regarded as a good quantum num- 
ber. The 77 = + spin ground state |+) then has J a = 0, 
whereas the state |— ) lives in a degenerate doublet space 
with J fj — ±2. This then raises the corresponding spin 
energy of the |— ) state by a term of order v a /L, as shown 
explicitly in Eq. (|35|l . As a consequence, the 77 = — 



branches are shifted upward with respect to the 77 = + 
branches (Fig. |2t>)) which destroys the hc/2e periodicity 
of£($). 

The spin current carrying i] — — states may (but need 
not) be shifted up in energy so much that $ = hc/2e 
ceases to be a metastable minimum of £ ($). This is just 
the case for a non-interacting system. In the case where 
the uplifting of the 77 = — state is not as large, $ = 
hc/2e persists to be a metastable minimum in the energy 
versus flux curve. The difference between the ground 
state energy at $ = and $ = hc/2e is thus given by: 



|£b(*o/2) - Eo(0)\ = - min ( v<r , K p v p ) 
(isotropic spin) 



(55) 



Interestingly, Eq. I|55|) provides information about the 
Luttinger parameters of the spin sector. It has long 
been known that for a Luttinger liquid, the Luttinger 
parameters can be determined from of the ground state 
properties^ This technique is often applied to infer the 
charge Luttinger parameters, e. g. by calculating the 
ground state energy as a function of particle density and 
magnetic fluxiSi Eq. H55|) shows that the same technique 
may be used to infer spin Luttinger parameters, provided 
that v a < K p v p holds. In SU{2) invariant systems, the 
spinon-velocity v a may thus be obtained. Note that in 
this case, the J a = ±2, N a — states corresponding to 
the 77 = — branches in Fig. |2Jd) are degenerate with states 
having J a = 0, N a = ±2, which carry no spin current 
but have a net azimuthal spin projection S z = ±1. This 
degeneracy follows from Eq. I|35ll with K a = 1. The lat- 
ter states, however, will generally be lower in energy for 
spin gapless systems without SU(2) invariance. This fol- 
lows because one has K a > 1 in this case, since K a < 1 
would always lead to a spin gap. We thus predict that 
the branches corresponding to the metastable minima in 
Fig. Et)), if present, will have a net spin, rather than a 
net spin current, in models without SU(2) invariance. In 
this case, v a in Eq. I|55l) is to be replaced by v a /K a w& 

When N = 4m, all the patterns in Fig. [21 are shifted 
horizontally by hc/2e. In this case the global minima in 
Fig. |2h) are located at odd multiples of hc/2e. Hence 
the function £ (<£>) can distinguish the cases N = 4m and 
N = 4m + 2 in the case of gapless spins, but not in the 
case of gapped spins. The same result had also been 
observed for the t-J-J' models 

We now turn to the case of odd particle number 
N = 2m + 1 . In this case the selection rules Eq. (|llafl re- 
quires both N a and J a to be odd, reflecting the fact that 
there must be a dangling spin. (Of course, with a dan- 
gling spin the system cannot have a spin gap.) Now the 
quantum number 77 = exp (— in J a /2) may take the val- 
ues ±7. The two corresponding subsets of the spin state 
space are related by the transformation J„ — > —J a , which 
leaves the Hamiltonian invariant. As a consequence, the 
spin ground states \n = ±i) are exactly degenerate, and 
an exact hc/2e periodicity is obtained for the Hamilto- 
nian (|12|) at any system size, regardless of whether Hi is 
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-0.5 0.5 g,yg, -0.5 0.5 



FIG. 2: Energy branches as function of flux for even total particle number N, with and without spin gap. £($) is given by 
the lower envelope. The alternating label r\ = ± describes the spin state corresponding to each branch, where N = 4m + 2 is 
assumed, a) Spin gapped case. The flux period is $o/2. b) No spin gap. The $o/2 flux periodicity is destroyed by a relative 
shift between the rj = + and r\ — — branches. 



relevant or not. Also, since J p is now odd as well, the 
pattern shown in figurc|2Ji) will be shifted horizontally by 
hc/4e. Hence the global minima of will be located 
at odd multiples of hc/ie in this case. 

C. Discussion of the results at T — 

The results presented in the preceding section are ex- 
act for the Hamiltonian Eq. (|12|) . which is believed to 
be the low energy effective theory for all one-dimensional 
systems with gapless, linearly dispersing charge degrees 
of freedom^ These results may thus be expected to be 
representative for this entire universality class. To rig- 
orously justify this point, the effects of higher order, 
less relevant operators should be included into the model 
studied above. We will not carry out such a detailed 
analysis here. Rather, we will point out some expected 
modifications due to less relevant operators, and argue 
for the robustness of the basic results derived above by 
comparing them to special examples of microscopic mod- 
els, where the features of £($) are known analytically or 
numerically. 

The flux period of the repulsive Hubbard model was 
studied in Ref. |2?]. These results agree well with our 
findings for the spin gapless case. In particular, for odd 
particle number N the global minima of £($) are at odd 
multiples of hc/4e, and the flux period is hc/2e. While it 
may seem surprising that the flux period does not distin- 
guish the spin gapless, odd particle number case from the 
spin gapped case (except for the position of the minima) , 
the microscopic origin of the hc/2e period is of a rather 
different nature in the two cases. A more subtle effect 
may demonstrate this: If one calculates the £ ($) of free 
electrons for odd TV, one indeed finds that hc/2e is the 
flux period. However one also finds that there exist cor- 
rections to the hc/2e period at order 1/L in £ (<&). These 
corrections are due to the band curvature neglected in 
the Hamiltonian Eq. (|12f> . Similar corrections to 
also exist at odd N for the t-J-J' model. They can be 
calculated using the method discussed in Ref. Such 
corrections in powers of 1/L, however, were not found in 



the t-J-J' model for the spin gapped case, where only 
exponentially small corrections were observed. We thus 
argue that corrections to the hc/2e flux period generally 
scale as 1/L in the odd N case, while they are exponen- 
tially small in the spin gapped case. The behavior in 
the latter case can be attributed to the fact that the spin 
gap generally causes an exponentially small sensitivity to 
boundary conditions in the spin sector. This point will 
be further clarified in Appendix lAl 

The analysis in the preceding section predicts the 
ground state to be unique on the branches of £(&) which 
contain the global minima. We expect this to be obeyed 
by general Hamiltonians. However, the four- fold degen- 
eracy which we found between the N a = 0, J„ = ±2 
and N a = ±2, J a = states on the metastable branches 
in the isotropic, even N, gapless spin case is an artifact 
of our restriction to the Luttinger Hamiltonian Eq. J5J • 
Rather, the true eigenstates are given by a triplet and a 
singlet to be formed from these four states, giving rise to 
a small splitting. However, except for this effect, the im- 
plied degeneracies at the crossings between the branches 
remain valid: Although the conservation of J p is approx- 
imate once higher order operators are allowed, a change 
of J p by 2 implies a change of momentum by 2k f. Hence 
the states at a branch crossing will not be mixed, and the 
cusps in £($) will remain sharp for general models^ Fi- 
nite size studies of the Hubbard modeffi show that the 
patterns displayed in Fig. [5] indeed emerge very clearly 
in numerical simulations carried out at moderate system 
size, both for the spin gapped (attractive) and gapless 
(repulsive) case. 

We thus conclude that all systems which can be char- 
acterized as Luther-Emery liquids have the hc/2e flux pe- 
riod. In particular, deviations from the patterns in Fig. 
|2]such as the appearance of additional minima at higher 
fractions of a flux quantum must be attributed to finite 
size effects. Such additional minima at <&o/n are known 
to occur in the large [/-limit of the Hubbard model^i, or 
the small J-limit of t- J- type models^, for fixed system 
size. The criterion for such finite size effects to disappear 
is that the amplitude of the oscillations AE ~ 1/L from 
Eq. (|53|l is small compared to any other energy scale of 
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the system. In the above cases, the relevant competing 
scale is J ~ t 2 /U. The associated crossover is clearly ob- 
served in Ref. where the t-J z model is studied: For 
t/L < J z , the model displays the spin gapped behavior 
shown in Fig. This is a consequence of the Ising spin 
gap of this model. The similar crossover for the repulsive 
Hubbard model is shown in Ref. l2?l where the pattern of 
Fig. |2p) emerges (with the necessary shift for N = 4m). 
We note that the appearance of local minima separated 
by half a flux quantum from the global minima in the 
repulsive Hubbard model is sometimes interpreted as a 
sign of pairing. 30 We stress, however, that this case does 
not meet the criterion of a $0/2 flux period as we define 
it, since a small splitting of order t 2 /U remains between 
the two types of minima of £($), which does not vanish 
as the system size is taken to infinity. 

D. Non-zero temperatures 

Finally, we briefly comment on the expected general- 
ization of our findings to finite temperature. The be- 
havior stated below can be verified straightforwardly at 
the special solvable Luther-Emery point of the model 
Eq. 1(12)1 (see Appendix 0. For T > 0, we consider 
the modulations of the free energy F(T, <E>) as a function 
of flux. The particle number is held fixed, i.e., the aver- 
ages are taken in the canonical ensemble. (If the particle 
number were allowed to fluctuate, the even/odd effects 
discussed above would considerably weaken the sensitiv- 
ity to flux.) 

Under these conditions, the observations made above 
for the ground state energy will carry over to the free 
energy as long as T < AE = 0(1/ L). However, the limit 
of Eq. is not to be taken here, because the ampli- 
tude of the free energy modulations is proportional to 
exp (-const T/AE) rather than AE when T > AE. In 
the spin gapped, even particle number case it remains 
true that terms violating the hc/2e periodicity are sup- 
pressed by a factor ~ exp (—const A CT L/w CT ), where A a 
is the spin gap. Comparing the two exponential fac- 
tors, these hc/2e violating terms will be negligible un- 
til T is of the order of the spin gap. At this tempera- 
ture, the amplitude of £ ($) is already exponentially sup- 
pressed, provided that the system is large enough such 
that AE <C A CT is satisfied. This again shows that the 
hc/2e period will be obeyed as long as AE is the smallest 
energy scale (other than temperature) of the system. 

We note that the cusps between the branches in Fig. 
121 will be smoothened by thermal fluctuations, giving rise 
to a finite negative curvature and paramagnetic^ effects. 

V. CONCLUSIONS 

In this paper, we demonstrated that the ground state 
energy of Luther-Emery liquids will generally exhibit an 
hc/2e flux period. While this statement holds in a strict 



sense in the limit of large system size, finite size devia- 
tions are expected to be exponentially suppressed in the 
system size. This result had been anticipated in an ear- 
lier work on a particular microscopic realization of the 
Luther-Emery liquid, the t- J- J'-modelA Here, we gen- 
eralized the result of Ref. by showing that the hc/2e 
flux period is implied by the widely accepted low en- 
ergy effective theory describing such a phase. As a re- 
sult, we clarify why the state of the spin sector impacts 
upon the flux period when it is commonly believed that 
in one dimension spin and charge decouple at low ener- 
gies. An important aspect of our findings is that in sys- 
tems with even particle number N, the hc/2e period is 
triggered by the spin gap (i.e. pairing) alone and is inde- 
pendent of whether the superconducting pair-pair corre- 
lations are the dominant long-distance/time correlation 
function. This may be of particular value for the correct 
interpretation of numerical work. In addition, we have 
also discussed the expected finite temperature general- 
ization of our findings. 

In Ref. 0we stressed the SU(2) invariance of the model 
discussed there. This requirement has been relaxed in 
the present discussion, where we did not enforce SU(2) 
invariance. Instead, only the weaker requirement of a 
conserved z-component of the spin S z = N a /2 was found 
necessary. In the anisotropic case, we must also require 
that the spin gapped ground state has S z = (see foot- 
note l26|) . which should be automatic in the isotropic case. 

Our findings underline the intuitive notion that every 
spin gapped system with linearly dispersing charge modes 
should share some features of a superconductor. 
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APPENDIX A: REFERMIONIZATION OF THE 
SPIN HAMILTONIAN FOR FINITE SYSTEM 

It is well known that the sine-Gordon model 

h„ =^ J dx:^K a (d x e a (x)) 2 + ^(dMx)) 2 y. 

J dx:cos(40 CT (x)): , (Al) 

can be mapped onto the massive Thirring modelAS 3 . This 
mapping permits a rather direct demonstration of the ex- 
ponentially small energy difference between the 77 = + 
and 77 = — states discussed in section llV Al In addition, 
when K a = 1/2, i.e. at the Luther-Emery point, the 
massive Thirring model reduces to a massive free fermion 
Hamiltonian which allows the exact calculation of vari- 
ous physical quantities. In particular finite temperature 
results can be obtained at the Luther-Emery point easily. 
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In the notation established in section lHIl the mapping 
onto the massive Thirring model can be performed by the 
introduction of the following spinless fermion operators: 



where A r 



ij-rx I \ 

e ^4 . e i(2r<M*)+M*)J . 

VL r ' ' (A2) 



Here, the symbol : : denotes a normal ordering conven- 
tion analogous to that defined in section ITTT1 but where 
(pr lS is replaced by the field 



0*0 = \ E S ( 3l PrA x ) ~ <P-r,s( x ) 



= J E S ( 3< Pr,s ~ <P-r,s) = 2r( t>o 



(A3) 



It is interesting to note that the operators ip r (x), in 
terms of which the spin Hamiltonian Eq. IjAlfl is best ana- 
lyzed, lead out of the physical Hubert space H. p hy S '- Apart 
from affecting the spin current quantum number J a , they 
also change the total number of net excited spins, N a , by 
1. In the physical Hilbert space, such a change must al- 
ways go along with a change of charge quantum numbers, 
which are not affected by ij) r {x). It is quite natural that 
the action of a single "fractionalized" operator such as 
ip r {x) will lead out of the space of physical states. At 
this point the larger space TLf ra c discussed in Section ITTT1 
becomes indispensable, as it allows us to define operators 
such as ip r (x) in the first place. It is clear, however, that 
these operators enter the Hamiltonian only in appropri- 
ate pairs, which leave the physical subspace invariant. 

Using standard methods reviewed in Ref. 17, it is 
straightforward to show that the field defined in Eq. (|A2|) 
satisfies the required anticommutation relations. The 
additional factor exp(iirrx / L) in Eq. (|A2|) will be com- 
mented on below. For now we note that it gives rise to the 
following boundary conditions for the spinless fermion 
fields: 



ip r {x + L) = e 



m(N„+J a /2-l) 



i>r{x) 



(A4) 



As is relevant to section IIVI in the following we will 
concentrate on the case N a = and J a even. Eq. (|A4|) 
then tells us that the sector r\ — — is represented by 
fcrmions obeying periodic boundary conditions, and the 
rj = + sector by fermions obeying antiperiodic boundary 
conditions. Note that without the additional twist in 
Eq. HA2|> . it would have been vice versa. 

Using the methods discussed in Ref. IrH one may now 




FIG. 3: Ground state of the fermion Hamiltonian Eq. I|A5(I 
for g = 0. The linear branches indicate the massless case 
where both iV CT and J a are good quantum numbers. Crosses 
indicate allowed momenta in the rj = + sector, boxes those in 
the rj = — sector. The encircled areas indicate the occupied 
states on each branch for the ground states in the N a = 0, 
Jo- = sector (crosses), as well as the N a — 0, J a = +2 sector 
(boxes). The latter state is degenerate with the N a = 0, 
J a — —2 ground state, where the k = mode is occupied on 
the r — — branch instead. For non-zero mass and N a = 0, a 
gap opens up and the upper band remains empty while the 
lower band is completely filled, where the allowed momenta 
depend on r\ as shown. 



show the equivalence of Eq. (|A1|) and 



E 



E 



dx (—irv) : tl)l(x)d x ip r (x) : 



dx m{i) r : ip].{x)ip- r {x) 



*E 



H„ 



o 

const 



dx : tpl(x)ip r (x) : : %}j^_ r (x)ip- r (x) : 

(A5) 



when 



,33 



|gl±l 

27ra 



1 

1 



-4K, 



(A6) 



Here, the symbol : : is as defined in Eq. ©, but the 
vacuum state is now the N a = J a = state that is 
annihilated by the field Eq. l|A3p . 

When K a < 1 the sine-Gordon model is massive, 
and the physics of Eq. ljA5(l is given by massive spin- 
less fermions. In particular, the fermionic interaction g 
vanishes for K a — 1/2, and the Hamiltonian Eq. (|A5|) 
becomes that of a massive free fermion model. This is 
the special point identified by Luther and Emery^ At 
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the Luther-Emery point the fermion dispersion relation 
is given by e(k) — y/ (vk) 2 + to 2 . In a range of K a val- 
ues around 1/2, the g-term only gives rise to quantitative 
corrections. 

We now return to the factor exp(iirrx/L) in Eq. (|A2|I 
(or the boundary conditions Eq. $Aty ) and show that 
with this factor (or the boundary condition specified in 
Eq. i|A4|l ) the boson and fermion theories are consis- 
tent. To illustrate that we compare the ground state 
degeneracies in the non-interacting massless case for the 
bosonic and fermionic theories, i. e., we let K a = 1/2 
and gi± = 0, which results in m = g = 0. 

In the absence of gi± both N a and J CT are good quan- 
tum numbers. Let us denote the ground state in the 
N a , J a sector by \N a , J a ). From Eq. (|35l) we recall that 
|0, 0) is the non-degenerate global ground state of the spin 
sector, whereas the states |0, 2), |0, —2) form a degenerate 
doublet. That this also holds in the fermionic represen- 
tations of the model is just achieved by the boundary 
conditions Eq. I|A4(1 (see Fig. [3|and the caption). Note 
that since the fermions Eq. l|A2fl are derived in terms of 
bosons, the fermion occupancies claimed below Fig. 
follow by derivation, not by definition. In fact one may 
easily evaluate 



(N„ J„\4(k)cr(k)\N„, J a ), 



(A7) 



where 



cl (k) = —= / dxe 
VI Jo 



ikx Tf 



2l ir , , 

k = — n + -(N <7 + JJ2-1) 
1j Li 



(A8) 



by plugging in Eq. I|A2|1 . and verify that the occupancies 
identified for the various states in Fig. [3] are correct^ 

The spin sector of the model is now represented in 
terms of fermions satisfying a conventional boundary 
condition given by Eq. I|A4|I . The difference between the 
ground state energy for r/ — + and r\ = — thus becomes 
the change in the fermion ground state energy induced 
by a change of the boundary condition, or equivalently 
the modulation of the fermion ground state energy 
caused by an AB flux. When the sine-Gordon model is 
massive, the fermions form an insulating state. Then, 
the sensitivity of their ground state energy to the bound- 
ary condition will vanish exponentially with the system 
size, as is well known from the general arguments given 
by Kohn£ and Thouless^. The greatest advantage of 
the refermionization occurs at the Luther-Emery point. 
For in that case, a non-trivial interacting bosonic theory 
is mapped onto a free fermion theory. In particular, 
at the Luther-Emery point g = one obtains from a 
direct calculation that the energy difference between 
the ground states for r\ = + and r\ = — vanishes as 
mexp(—mL/v), as we claimed earlier. 



APPENDIX B: NUMBER AND PHASE 
VARIABLES FOR CONTINUUM AND LATTICE 
HILBERT SPACES 

The Hilbert space of the Luttinger Hamiltonian 
Eq. denoted as the "physical" Hilbert space Tt p h y s 
in the bulk of the paper, can be decomposed as 



H 



phys 



(Bl) 



Here, the spaces Ti hr s contain all the degrees of freedom 
associated with the bosonic excitation spectrum, whereas 
Ti.f Ir,s contains the degrees of freedom of the operator 
N r>s . The subscript "Z" stands for "lattice" and reminds 
us of the discrete nature of N ryS in the physical Hilbert 

space: The Hilbert space basis of Tif r,s is given by a set 
of non-degenerate eigenstates of N r>S) whose spectrum 
consists of all integer numbers. 

In the process of bosonization, however, we introduce 
new linear combinations A^y, J„ of the JV rja (Eq. ©). The 
spectrum of these new operators is likewise integer, yet 
not all possible combinations of integer eigenvalues are 
physically allowed. This "residual coupling" is not evi- 
dent from commutation relations, since the operators N Ul 
J„ all commute as the N r ^ s do. Hence, once we bosonizc 
an enlarged Hilbert space becomes more natural, where 
the eigenvalues of the operators N v , J„ are independent. 
This is the Hilbert space Hfrac of "fractional" excita- 
tions. The physical subspacc TL p h ys is then characterized 
by the fact that the selection rules Eq. I jllj l are satisfied. 
It is clear that in Tif rac , the spectrum of the N r<a must 
also contain certain fractional values. Formally, we find 
it convenient to introduce an even larger Hilbert space 
7i, where the spectrum of the N r ^ s is continuous. The 
benefit of this is that the conjugate phase Tp r s of these 
operators then becomes meaningful. This, in turn, al- 
lows us to construct unitary "ladder operators" which 
change the eigenvalue of the N rtS by arbitrary amounts. 
This formalism is of particular advantage in Appendix 
lAl where fractionalized spin fermion operators are con- 
structed. Below we present some fine details of this em- 
bedding of Hphys into the larger space H. 

For this purpose let us consider a single operator N 
and its conjugate variable (p such that 



if, N =i 



(B2) 



holds. An analogy to the quantum mechanics of a point 
particle moving in one dimension is obtained if we iden- 
tify N = x and (p = — p, where x and p are the coordinate 
and the momentum of the particle. In this context, it is 
familiar how to construct a Hilbert space TL N such that 
TV and ip are well defined on a dense set, Eq. (|B2|I is 
satisfied and the spectrum of N is unbounded: It is the 
Hilbert space of square integrable functions of the vari- 
able N. From the commutation relation Eq. I|B2(1 . it is 



12 



clear that the spectrum of both N and tp has to be con- 
tinuous and unbounded. In particular, we can construct 
shift operators exp(i<pa) satisfying 



N. 



i(pa 



a e 



iipa 



(B3) 



which shift the value of N by an arbitrary amount a. 

We note that, as is familiar from the point particle 
analogy, the "position" and "momentum" eigenkets | N) 
and | tp) are not strictly contained in the Hilbert space 
H N of "proper" vectors, but are "generalized" states in 
the usual sense. Here, we will not attempt to introduce a 
different notation for proper and generalized states, nor 
for the proper Hilbert space and its extension containing 
generalized states. We refer the reader to Ref. |3^ and 
references therein for details, and simply note that the 
kets | N) and | 0) satisfy 

(N' | N) = 6{N' - N) , (tp' | ip) = 2ir5{tp' - ip) 

1 

2^ 

(B4) 



\tp) = dNe 



-iipN 



N) AN) 



dtpe ivN \tp) 



Suppose now that a physical problem is defined on a 
subspace Ttf , which is given by the discrete "lattice" 
represented by the eigenkets | N) for integer N. Since 
these kets from a countable Hilbert space basis in H(* , 
it is natural and convenient to introduce a new scalar 
product on Tif via: 



,{N' | N) t = 5 N >,N 



(B5) 



This differs from the scalar product in TL N only by an 
infinite multiplicative factor. In the above, | N), denotes 
the same vector as | N) , but endowed with a different 
scalar product. Eq. (|B5|) means that the | N) for integer 
iV become a complete orthogonal set of proper vectors 
within TiV . Within %F , one may now define a "crystal 
momentum" operator tpi, whose eigenkets are defined to 
be 



|V>,= E 



-itpN i 



(B6) 



One observes that these eigenkets are periodic in tp with 
period 2n, hence for dcfmitcncss the eigenvalues must 
be restricted to lie within the "Brillouin zone" (— tt, 7r], 
where 



l {<p'\(p), = 2ir6((p' -(f), tp', tp e (— tt, 7r] 



(B7) 



holds. If we now denote the restriction of N to Tif by 
Ni, we find that the commutator [tpi, N[] is not quite anal- 
ogous to Eq. i|B2(l . This has been examined in detail 
by Sch6nhammer^2i2& However, for the applications we 
have in mind here, this difference never matters. This is 
so since all physical observables, including the Hamilto- 
nian, depend on tpi only via integer powers of exp(iipi), 
and since the equations 



«v 



N) = \N+ 1) 



\N)i = \N + l), 



(B8) 



hold. Hence exp(itpi) and exp(i^>) act identically on Ti.f , 
and we may express all observables in terms of either of 
these operators. By means of Eq. (|B6|> . the kets | tp), £ 



(B9) 



TC^ are identified with the following kets in Tt N : 



tp\ = J2e-^ N \N) = J2\v> + 2nn) 



Note that the norm of the right hand side with the scalar 
product in 7i N formally computes tooox (5(0), which by 
comparison with Eq. I|B7|I is larger by an infinite mul- 
tiplicative constant than h {tp \ tp) v Recall that the same 
relation also holds for (N | N) and t (N | N) l by definition 
of the scalar product in Ttf . Eq. (|B9J| makes it clear 
that the kets | tp), G Tif are not to be identified with the 
kets | tp) E 7i . The latter are not periodic in tp, and 
cannot be constructed solely from kets | N) with integer 
N. However, 



{tp + 2vm | N) = t (tp | N) t 



Vn e Z 



(BIO) 



holds. Thus whenever it is clear from the context that we 
are working in 7i¥ , we may drop all labels keeping 
the periodicity of the states | tp) £ Tif in mind. 

We can now define the Hilbert space Ti introduced in 
the main part of the paper as 



H = H Nr ' B ® U h 



(Bll) 



The advantage of embedding TL p hys into the space Ti is 
that now the space Hf rac , satisfying 7i p hy S C Hfrac C 7i, 
can be generated easily through the action of the oper- 
ators exp(j#„), exp(i<^„) on TL p h ys . The bookkeeping is 
greatly simplified by the simple commutation relation of 
Tp r s and N r ^ s , Eq. fflljl. valid on TL. We finally note that 
Tt f r ac can be written as a product analogous to Eq. (|B1|I , 
involving a space containing bosonic degrees of freedom 
and four discrete "lattice" spaces containing the degrees 
of freedom of the quantum numbers N v , J„. All the 
above therefore holds in an analogous way for TLfrac, and 
in particular lattice versions of the phase operators <^„, 
9 V can be constructed if desired. 

In Ref. |H, it has been noted that the use of the canon- 
ical commutation relations Eqs. 1)20(1 . I(B2(I in construct- 
ing the Klein factors does indeed yield correct results 
in bosonization, even though these relations cannot be 
rigorously justified for operators that are restricted to a 
discrete space. We believe that the embedding procedure 
discussed here provides a proper explanation for this ob- 
servation. 
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